35  603  THE  INUERCE  BACK  CCATTERINC  PROBLEM  FOR  THE 

iCHROEHNCER  EQUATION  IN  TWO., 'Ll;1  NEW  YORK  LINIU  NY 
COUFANT  INST  OF  MATHEMATICAL  CCIENCEC  C  C  MO  Pa  WET  2 

:::fied  01  app  3?  N00014-35-K-0195  f  :  12/1 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-  1%3-A 


AD-A185  603  > 

■  ■ 

% 


2b.  OECLASSIFICATION/OOWNG 


4.  PERFORMING  ORGANIZATION  HfcPORT  NUM8ER(S) 

o&D- 


6.  NAME  OF  PERFORMING  ORGANIZATION 

COURANT  INSTITUTE  OF 
MATHEMATICAL  SCIENCES 


6c.  AOORESS  IClly.  Slate  and  ZIP  Code) 

NEW  YORK  UNIVERSITY 
251  MERCER  STREET 
NEW  YORK,  NEW  YO 


8*.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 


8c.  AOORESS  idly.  Stale  and  ZIP  Code) 


•  ^  • 


^^DOCUMENTATION  PAGE 

lib.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION/ 


ON  STATEMENT  A 


Approved  for  public  releasel 
Distribution  Unlimited 


7a.  NAME  OF  MONITORING  ORGANIZATION 

DEPARTMENT  OF  THE  NAVY 
OFFICE  OF  NAVAL  RESEARCH 


7b.  AOORESS  (City,  State  ond  ‘/IP  Code) 

800  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable  ^JQOO  - 

N-Oflfl  1 f-8  5 -K-0 1 9  5 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 


1 1.  TITLE  (Include  Security  Classification)  InVfiTSe  BflCk 

Problem  for  the  Schroedinger 


PROJECT 

TASK 

WORK  UNIT 

NO. 

NO. 

NO. 

ithleen  S. 

Morawetz 

13«,  TYPE  OF  REPORT 

Interim 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM _  TO 


14.  DATE  OF  REPORT  IYr..  Mo..  Day )  IS.  PAGE  COUNT 

April  1,  1987  26 


COMMUNICATIONS  IN  PURE  AND  APPLIED  MATHEMATICS 


18.  SUBJECT  TERMS  { Continue  on  reverse  if  necessary  and  identify  by  block  number) 


17. 

COSATI  COOES  1 

field  ! 

GROUP 

SUB.  GR.  1 

Inverse  Scattering 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  bloc*  numberT^*-,.  ,  ,  ,  j  .  , 

This  paper  is  concerned  with  determining 

a  potential  q(x)  for  the  steady  state  Schroedinger  equation  in  two  space  variables, 
x  -  (x1,x2): 

1..1  (A  +  U)2  -q)u  =  0  . 

It  is  assumed  there  are  no  bound  states. 

The  data  we  are  given  come  from  the  far  field  scattered  by  plane  waves  impinging  in 
a  range  of  directions  but  measured  only  in  the  opposite  directions.  More  precisely,  let 
p(e,03,x)  be  a  solution  of  (1.1)  which  for  e#x  «  -  00  behaves  like 

eiu)e  x^i  ^  Here  e  is  a  unit  vector. 


IBUTION/AVAILABILITY  OF  ABSTRACT 
UNCLASSIFIEO/UNUMITEO  CD  SAME  AS  RPT.  □  OTIC  USERS  □ 


22s.  NAME  OF  RESPONSIBLE  INDIVIDUAL 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22b.  TELEPHONE  NUMBER  22c.  OFFICE  SYMBOL 

( Include  Area  Code ) 


1-212-998-3200 


EDITION  OP  1  JAN  73  It  OBSOLETE. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


**  5  20 


1  23 


JCathleenS^Morawetz 
DO  FORM  1473, 83  APR 


We  consider  the  space  (£  of  bounded  functions  also  lying  in  L*  and  assume  q 
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where  s^  is  called  the  far  field.  The  back  scattered  field  is  measured  in  the 

direction  -e  and 

1*3  s(o),e)  *  s^e.u^-e) 

The  problem  under  study  is  the  determination  of  q  from  s(co,e).  A  theorem  is 
proved  here  that  if  s  is  sufficiently  small  so  that  the  Born  approximation  is 
small  in  norm  in  the  space  (£  there  exists  a  unique  solution  q  to  the  back 
scattering  problem. 


The  Inverse  Back  Scattering  Problem  for  the  Schroedinger  Equation 

in  Two  Space  Dimensions 
Cathleen  S.  Morawetz 


1  lntXfldUfitiflD 

This  paper  is  concerned  with  determining  a  potential  q(x)  for  the 

steady  state  Schroedinger  equation  in  two  space  variables,  x  -  (x^,x2): 

(1.1)  (A  +  w2  -  q)u  -  0  . 

lc  is  assusMd  chars  srs  no  bound  scscss.  t- 

The  data  we-aee  givey  come  from  the  far  field  scattered  by  plane  waves 
impinging  in  a  range  of  directions  but  measured  only  in  the  oppo.^1  re¬ 
directions.  More  precisely  let  p(e»u,x)  be  a  solution  of  (1.1)  which  tor 
csX  -  -  «  behaves  like  m^9****  plus  a  scattered  wave  behaving  like  e*we#A 

Here  e  is  a  unit  vector.  Jlearly  this  requires  some  conditions  of  decay  on 

q- 

We  consider  the  space  Q  of  bounded  functions  also  lying  in  L1  and 
assume  q  «  Q.  This  is  a  less  than  perfect  space  from  the  point  of  view  ot 

phvsics  since  it  does  not  allow  for  singularities  in  the  potential.  Howewi 

in  other  work  [1]  a  study  is  made  of  an  iteration  for  finding  a  smow 

potential  from  its  back  scattered  field  and  this  space  suffices  for 

study,  see  Appendix  2.  One  would  prefer  for  mjuserical  purposes  to  conswu 
only  compactly  supported  potentials  but  the  theory  does  not  work  for  thi  ^ 
It  should  be  noted  that  the  2D  problem  is  barely  admissible  for  computation 
as  it  Involves  four  independent  varlbles  (e,w,x)  and  an  iteration. 

The  solution  p  of  (1.1)  is  written  as  e  ♦  p#  where  for  |x| 
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where  is  called  the  far  field.  The  back  scattered  field  is  measured  in 

the  direction  -e  and 

(13)  s(«>,e)  -  s^e.w.-e)  . 

The  problem  under  study  is  the  determination  of  q  from  s(u>,e). 

From  the  physicist's  point  of  view  this  may  seem  an  unreasonable 
restriction  on  the  measurement.  At  the  very  least  one  could  measure  the 
scattered  field  over  a  range  of  angles.  However  this  would  give  us  an  over- 
determined  and  hence  Intrinsically  unstable  problem.  Ue  would  be  giving  a 
function  of  two  angles  and  the  frequency  u>  which  is  too  much  data  Iol 
determining  q.  a  function  in  polar  coordinates  of  one  scalar  and  one  angle 
It  is  reasonable  to  expect  that  the  linearized  theory  (the  Born  approxim a> 
tion)  determines  how  much  data  should  be  prescribed,  tfe  check  this  now  in 
2D.  First,  rederiving  the  Born  approximation  and  then  proceeding  to  the 
main  result: 

For  sufficiently  small  values  of  s(e.ui),  there  is  a  unique  potential 
if  there  are  no  bound  states. 

By  the  obvious  linearization,  if  is  the  Born  approximation  to 
scattered  field  p  then 

/  m.  .  2X  -iwe*x 

(A  ♦  w  )pR  -  q  e 

Hence  introducing  the  fundamental  outgoing  solution,  a  Hankel  function, 
we  have  up  to  a  constant 

PB(x)  -/  H(w|x*x'|)  q(x')  e“lu*#X  |dx'| 

For  large  distances  e*w*x^  |x|  *^PB(x)  tends  in  the  direction  -e  to 
J  v/l/"2  exp  iw( -e|x|+x' )  •  e*x'  -  |x|)  q(x')  dx' 

-  «*l/2  J  q(x<)  dx'  , 

up  to  a  ph««t  factor.  Thus  totting  q  -  q^  to  donoto  tho  Eorn  approximation 


( 1  4a) 


•  (#,<•»)  -  o  u 


qB(-2wo) 


where  is  ths  Fourier  transform  of  q  and  c  is  soms  numerical  factor.  Thus 
the  relation  betveen  the  approximate  potential  and  the  back  scattered  fieLd 
can  be  read  off  by  the  Fourier  transform.  If  q^  lies  in  Q  we  note  that  s  is 
bounded  but  more  than  that  is  required  of  s  to  determine  q.  It  is 
reasonable  to  suppose  that  if  the  scattered  field  is  small  in  a  suitable 
norm  such  as  that  qR  is  small  then  the  exact  q  can  be  found  as  a  convergent 
series  in  powers  of  qR.  This  has  been  done  for  3D  by  Prosser  in  a  series  of 
papers  [2].  See  also  Newton  [3]  where  the  full  scattered  field  is  involved 
and  Fadeev  [4] . 

The  use  of  the  Radon  transform  approach  in  3D  is  analyzed  by  Yagle  et 
al.  [5].  Here  we  have  transformed  away  from  the  Radon  transform  which  is 
used  in  [ 1 ] . 

The  approach  used  here  utilizes  the  underlying  formulation  of  [1]  and 
regards  the  given  data  as  q^(x) .  Then 

a.-*)  q-(x)  -  q(x)  ♦  c  /  |u>|  p  (e.w.x')  q(x' )  el{^*(2x'x  >  |  dx '  |  dw  d* 

o  s 

where  c  here  and  elsewhere  is  a  generic  real  constant  and  - 
cos  $,  sin  0).  The  derivation  is  in  appendix  1. 

The  formal  expansion  is 

m 

(1.5)  q(x)  -  £  q  (x)  , 

1  n 

where 

(16)  qt(x)  -  qR(x) 

q  (x)  -  c  /  M  p  (o.w.x*)  q(x')  *  du  d$  , 

n  n 

jnd 

(17)  APn  ♦  *  qpn-l 

•  iu>e*x 

P0  -  • 

Ua  would  Ilk*  to  prova  that  (1.5  -  1.7)  load  to  a  unlqua  aolutlon  toi  , 
<  Q  In  taraa  of  q^  <  Q. 


Certain  aspects  of  this  problem  are  more  difficult  than  in  the  ilJ 
problem.  This  is  most  easily  seen  from  the  corresponding  wave  equation 
Denoting  the  corresponding  variables  by  capital  letters,  let 


(1.8) 

P(«,t,x')  -  c  /  p(«tW>x')  d Im 

Then 

(1.1*) 

Ptt  '  AP  +  qP  “  0  ’ 

P  -  0  for  t  -  •  •  x  <  0  , 

and 

P  “  S( t-«»x)  a»  t  -»  -  •  . 

Then 

,  _  ★ 

ll.2) 

Pstt  -  AP^  ♦  qP^  -  -  q4 (t-«*x) 

(1.3)* 

S(«.t)  -  ic  /  «lwt  w'l/2  qB(2w«)  du 

(1.4)* 

Ait  d? . 

qB<x>  ~  q(x)  ♦  c  /  jjf  (•,««(2x 

-  a,  x1)  q(x')  |dx' | 

£P 

06)  qn  -  c  J  ^ ^  ^  («,«• (2 x-x* ) -a ,x* )  q(x')  |dx'| 

,  l  7)#  Pntt  *  APn  *  -  qPn-l 

-uh  po  •  <<«-••*>. 

The  dimension  problem  shows  up  in  (1.7)  .  The  nonhomo geneous  term  tu* 
n  >  2  is  not  of  compact  support  in  time.  A  crude  estimate  suggests  that 
there  is  logarithmic  growth  in  time  for  Pn  if  there  are  two  dimensions  in 
space.  In  three  dimensions,  by  contrast,  Huyghen's  principle  carries  the 
energy  of  P  rapidly  out  of  the  main  support  of  q  and  estimates  are  general iv 
easier  to  make. 

We  use  only  some  the  properties  of  q^  that  can  be  derived  troir 
■ i ire  -  formulation  (*)  and  then  use  direct  frequency  domain  estimates 
de fine 

a. 9)  I  I q I  I  -  lq|.  ♦  /  |q|  l<lx| 
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and  show 

(110)  ||qnll  s  cn  | )q| |°  . 

Thus  the  series  (1.5)  converges  for  ||q||  less  than  some  constant  and  can  be 
uniquely  inverted  if  | Iq^l |  is  sufficiently  small.  Thus  we  will  have  proved 
Theorem.  If  the  norm  of  the  Bom  approximate  ||q^||  is  sufficiently 

small  there  is  a  unique  potential  with  the  same  back  scattered  field  as  the 
Born  approximate  has  by  linear  theory. 

The  proof  of  (1.10)  consists  of  technical  (mainly  calculus)  estimates 
using  iterated  kernels  in  both  the  time  and  frequency  form.  In  section  2  we 
examine  which  turns  out  to  have  the  interesting  form 
(1.11)  J  d *'  Q(«*x,e)  Q(®' *x,e' )  «(••«')  , 

where  &  Is  the  Dirac  delta  function  and  Q  is  the  Radon  transform  of  q  at  x 

in  the  direction  perpendicular  to  e  -  (cos  $t  sin  9 ).  We  also  show  that  it 

2 

is  bounded,  |  1 I  I  *  c  I  I q I  I  • 

In  section  3  we  show  that  for  n  >  4,  q  is  correctly  bounded  in  L  :  in 

n  « 

section  4  that  q^  and  q^  are  also  correctly  bounded  In  L^.  In  section  :> .  ..v 
show  that  q^  lies  in  L'  for  n  £  3  again  with  the  right  bound.  The  proof  ot 
the  theorem  is  then  easily  completed. 

We  shall  ignore  constant  factors  to  the  power  n  as  they  influence  only 
the  radius  of  convergence  of  the  series.  However  we  shall  point  them  out 
where  they  occur. 

In  Appendix  1  we  rederive  the  basic  formula  (1.4)  and  (1.4)*.  In 
Appendix  2  we  show  how  an  iteration  to  achieve  the  solution  converges  usLi:>; 
previously  obtained  estimates. 
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2.  Th?  sags  n  -  2. 

To  derive  the  formula  (1.12)  we  may  assume  q  is  smooth. 


We  have  from  (1.7) 


?1  -  -  S  H(w|x-x'|)  q  e’iwe#x  |dx'|  . 


Thus  we  find  from  (1.6); 


where 


q o  "  /  K(x,x',x")  q(x' )  q(x")  |dx' |  |dx"| 


K  -  Re  c  /  |w|  H(w|x'-x"|)  elwc' (2x-x' -x") 


r  iw  (2x~x'  -x")  ..  T  ,  M x 

J  e  d*  -  c  J0(w|2x-x#-xM|)  , 

where  Jq  is  the  regular  Bessel  function  of  zero  order.  Hence  we  write  K  as 
K(pQf/>l)  where  pQ  -  |2x-x'«x"|,  p^  -  |x'-x"|  and 

K  -  Re  c  J  |w|  JQ(wpQ)  H(wp^)  do)  . 

-  eo 

This  integral  has  to  be  properly  interpreted.  First  we  have 

00 

K  -  Re  c  /  w  J0(w/>0)  OKw^)  +  HC-wp^))  dw  . 

To  be  outgoing,  H  is  the  Fourier  transform  of  a  function  of  t-|x|  and  hence 
behaves  like  e  at  •,  It  also  has  a  real  logarithmic  singularity. 

Hence 


(2.1) 


H(f )  -  i  Jn(f)  +  Nn(f)  . 


K  -  c  f  v  JQ(u >pQ)  ^(wp^  d<i>  . 

2 

From  Bessel's  equation,  (C-Jq)  +  {"Jq  -  0  we  find  +  Pq1^  “  0  and  hence 

wJN  -  (wJ^N  -  wJNw)w  with  f  -  Jq(wPq)  ,  N  -  .  Hence 

k  -  -f-j  (li»  -  ptjft)  +  J)  . 

P0-Pi  0— 
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We  have  used  the  logarithmic  behavior  of  N  at  w  -  0  and  the  regularity  of  J. 
(See  e.g.  Courant-Hilbert  Vol.  1,  p.  501  for  the  constants.) 

We  now  use  the  asymptotic  behavior  of  the  Bessel  functions: 


fi(P0>JN  -  Pjjfi)  -  J  (-  J  J*-  (cos  o (.Pq-Pi)  -  cos  O(p0-p1  -  |)) 


-  J  ~r  (cos  0(Po'pl^  +  cos  n(-Pn'P-\  ' 


0  ' 1  2' 


Hence 


K  -  lim  ~2~2  ‘  cos  n(^o'pl>>  +  O  *  2  7  ~p^  '2  7  T^)  c°s  CKPq-p i> 

£>«o  pQ-p^  1  0 

-  (  J  ~~  -  J  )  sin  0(p0+p.)) 

P\  p0  01 

which  we  write  as  lim  (K-+K-+K-). 

It  is  easy  to  see  that  the  last  term  behaves  in  the  limit  like  a  <5 - 
function,  in  fact 

K3  -»  c^Pq+p^)  l(p0Pi)  1/12  st.pQ+p^  . 

We  recollect  that  both  Pq  and  p^  are  nonnegative  when  this  factor  in 
the  kernel  is  applied  to  q(x’)q(x")  |dx’|  |dx"|.  We  can  show  there  is  no 

contribution.  First,  set  x  -  0  for  convenience.  Then  2x'  -  p^e  +  p^e',  2x" 
-  pQe  -  pj^e'  with  e  -  (cos  8,  sin  8),  e'  -  (cos  8’,  sin  8’).  The  integral 
becomes 

/  «(P0+P1)(P0+P1)‘1(P0P1)1/2  q(P0e  +  PXC')  q(P0e  -  P^')  dpQ  dpx  d 8  d 8' 

Set  pQ  -  R  cos  <t>,  px  -  R  sin  <t>  with  0  s  <f>  <  J  and  the  integral  is  the 
principal  value  of 

c  /  5 (R  cos  (^  *  £))  (cos(*-  2"))’1  sin1^2^  cosl^2<f>  q(R(e  cos  <t>  +  e'sin  <j>)) 
q(R(e  cos  ^  -  e'sin  ^))  R  dR  cty  d&  dd' 
which  is  0  since  cos  (<£  -  *)  *0. 
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For  K2  we  make  the  same  change  of  variable  and  we  have 
1/2  1/2  1/2  1/2 

RK2  dR  4  -  SM7//'stn1//  ^  *rsaa.L*  cos(0R  cos(,+  R  ^  . 

cos7^+sln/z^  cos -  *)  4 

For  R  cos(^  +  J)  >  5  the  integral  with  this  measure  over  q(x')q(x")  will 

tend  to  zero  as  fl  »  by  the  Riemann  Lebesgue  lemma.  For  R  cos(<£+  j)  <  <5  we 

note  that  R^dR  d^  is  antisymmetric  and  vanishes  for  5  -  0  so  again  the 

contribution  to  the  integral  can  be  shown  to  be  arbitrarily  small. 

Finally  we  turn  to  .  The  contribution  to  the  integral  is  again  found 

with  the  variables  pQ,  p^f  0,  0#.  Then  we  set  -  *?.  PQ-P  1  -  f  and  the 

integral  over  q(x')q(x")  |dx'|  |dx"|  becomes 

X  (1  -  cos  0»?)  q(^  r?(e+e' )  +  J  £(e-e'))  q(^  f?(e-e’) 

-  J  e(e+?'))(£2-*?2)  d£  d,?  d9  d*' 

or  with  Clrj  -  r ,  the  integral  is 

X  l  f  dT  X  |  ^(“4^  («+«')  +  4  £<•-«')  q(-£jj  (e-e') 

1  2 

-  T  £ (e+e ' ) ) (£  -  *r)  d£  d0  d0' 

n 

On  letting  0  *+  «>  we  have 

X  ^99  r  dr  f  (  q(J  £(fe-e' ))  q(-  J  ?(e+e'))  d£  dff  d* ' 

The  first  integral  may  be  replaced  by  c  and  we  are  left  with 
c  f  £  q(J  £ (e-e ' ) )  q(-  J  £(e+e'))  d£  d0  d0' 

But  (e-e' ) • (e+ef )  -  0  so  we  may  set  ^  f(e-e')  -  rE  with  E  -  (cos  <j> ,  sin  <f>) 
and  -  ^  £(e+e')  -  sE  where  E^E*  -  0  and  r^+s^  "  g  ^  so  £  d£  d0  d0 ' 
-  dr  ds  d4  d0  . 

We  may  then  write  the  contribution  to  q2  from  as 

-  t 

c  J  q(rE)  q(sE  )  dr  ds  cty  , 

which  is  the  same  as  (1.11).  This  form  will  not  change  if  we  let  q  belong 
to  the  wider  space  Q  provided  the  integration  can  be  performed. 


a 


! 


f 
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To  consolidate  we  write  this  result  as  part  (a)  of  the  following  lemma 
(shifting  x  back  from  the  origin). 

Lemma  2.1.  If  q  «  Q  then 

(a)  q 2  -  c  J  q(x  +  re)  q(x  +  se* )  dr  ds  d0 ,  and 

(b)  ||q2l|  <  c | | q | | 2 . 

Proof  of  (b) :  First  for  the  neighborhood  of  r  +s^  we  bound  by  the  L 
norm.  Then  we  set  r  -  R  cos  s  -  R  sin  0  and  x  -  0;  then 

q2  - 

c/  q(R  cos  <£  cos  0,R  cos  sin  0)q(R  sin  ^  sin  0,-R  sin  ^  cos  6)R  dR  d<t>  dv 
so 

Kol^clql  /  q(R  cos  4  cos  8%  R  cos  4  sin  8) 

[ cos  <f>\>6 

.  i  m  £<>§.  £l  ^  de 

COS  ^ 

+  cl^loo  J*  q(-R  sin  <j>  sin  9,  R  sin  4>  cos  Q) 

| sin  ^|>6 

.  R  ala  k  d(R  Sfn  4)  ^  d6 
sin2^ 

S  clql^Hql^/62  S  c |  | q |  | 2  . 

Secondly  for  the  l  estimate  we  calculate 

f  | q |  |dx|  £  c  /  /  |q(x+re)|  |q(x+se  )|  dx  ds  d8  d<£  . 

We  let  9  range  from  0  to  ir  and  r  and  s  from  -  ®  to  +  «  .  Set  x  -  (e*x)e  + 
(e^-X)e^  and  we  have  (for  convenience  dropping  the  absolute  value  signs), 
f  q((r+e«x)e  +  (e^»x)e^)  q((e*x)e  +  (e,*x+s)et)  dx^dx^dr  ds  d 8  . 

Taking  X£  -  e,*x,  x^  -  e*x  we  have  to  integrate  with  respect  to  8  the 

integral 


so  we  have 


/  q(r+xltx2)  q(x1>X2+s)  dx^dx2<lr  ds 

But  X  q(xljX2+s)  dXjds  -  I  |  q  |  |  ^  and  X  q(r+x1(x2)  dj^dr  -  ||qj|1 

X  |qnl  |dx|  s  c||q||J  . 

For  use  In  Appendix  2,  we  note,  as  is  easily  proved, 

Lemma  2.2.  If  q^  €  Q,  then 

5  "  X  qx(x  +  re)  q2(x  +  re')  dr  ds  d 6 
lies  in  Q  and  ||q||  <  c||ql||  ||q  ||. 
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3.  bounds  for  qn>  n  >  4. 

We  first  derive  some  estimates  above  Hankel  kernels  that  we  use 
frequently. 

Lemma  3.1. 

(a)  |  f  H(upi)  H(w^i+1)  qCx1)  Idx1!  |  <  c| |q| | ( | log  w|2  +  1), 
and  for  |u>|  >  0, 

(b)  |  f  H(WPi)  H(Wpi+1)  q(xL)  Idx1!  |  <  ||q||  . 

u  ,  i  i+1 . 

here  p^  -  |x  -x  | . 

Proof:  From  (2.1)  using  the  behavior  of  the  Hankel  function  we  have 

for  |  co  |  and  ^ess  t^an  k  say, 

|H(wpi)  H(w/>i+1)|  <  c  ( |  log  <p|  +  |  log  ^IKIlog  w|  +  |  log  Pi+1l)  • 

Hence  the  integral  of  (a)  over  |w|p^,  |a>|p^+^  <  k  is  less  than 
c  J  (llog  U>|2  +  |log  | 2  +  |log  Pt+1|2)  q(xl)  Idx1! 

<  c  ( |  |q|  |  x  |log  o,|  2  +  ||  q  |  |  L  |  log  5  I2  +  llqlljlog  5|2  S2) 

<  C  ||  q  |  |  (|  log  o,|2  +  1) 

For  |  cu  |  p  ^  >  k,  lwlPj_+^  <  k  we  obtain  a  similar  estimate  since  the  first 
Hankel  function  is  bounded.  For  |wp^|,  lwPi+]J  >  ^  we  get  the  bound  ||q||^. 
This  completes  the  proof  of  (a) . 

-1/2 

To  prove  (b)  we  note  that  for  |u>p^|  >  k,  |H(wp^)|  <  c  | cup ^  '  .  We 

have 

2  |  J  H(o,p.)H(o,p.+1)  quS  Idx1 1  |  <  J  (  |H2(o,pi)  |  2+  |  H(o,p  .+1 )  |  2)  |  q(x' )  |  dx ' 

kl*il  1  2 

We  estimate  the  first  term  by  o |  | q |  |  J  |H(wp.)|  P-  dp  + 

00  0  111 

clo,!'1  f  -1-  |q(xi)|  Idxl1  <  c|  |q|  |J<pf2  +  C|u>|  '1|  |q|  I  <  |  |  ~  1 1  |  q  |  |  Here 

^i 

we  used  |  f  p*1  q  dx|  <  R| |q| +  R  ^l  |q|  1^  for  any  R.  The  second  term  is 
essentially  the  same.  Thus  for  some  Q,  u>  >  ft  the  inequality  (b)  will  hold. 
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This  lemma  implies  that  we  can  obtain  estimates  for  integrals  with  many 

products.  Thus,  to  be  sufficiently  general  for  later  purposes 

-  2  - 1 

Lemma  3.2.  Let  H(u>)  -  min(|log  a>|  +  1,  |w|  )  and  J(~)  = 

min(l,  M'1/2).  Let  q  -  /  J(«p  )  II  H(wp)  II  qi(xi)  |dxl|.  Then  |q  |  < 

n  u  i-1  1  i-1 

cj(w)  ii(w)n/2  n  I  I  q1 1  I  . 
i-1 

Proof:  It  follows  from  Lemma  3.1  that  for  the  variables  xm,  m  *  1  and 

n  that  occur  in  pairs  in  the  product  of  the  Hankel  functions  in  the  integral 
one  can  estimate;  m  -  1,  n  are  exceptions  as  they  occur  in  the  J (o>Pq)  and  in 
H(u;pn)  .  There  are  n-1  good  terms  and  that  gives  (n-l)/2  if  n  is  odd  and  (n- 
2)/2  if  n  is  even.  This  leaves  "left  over  q*""  terms  but  there  is  no 
dependent  factor  for  n  odd  except  those  that  occur  in  the  argument  u>p^  of 
the  Bessel  function  J,  namely  x#  and  xn.  Thus  the  terms  in  can  be 
estimated  by  L^  integrals  except  for  two.  The  designated  integral  is 
bounded  by 


(3.1)  (H(u>) ) (n"1)/2  n  | | q1 | |  S  |J(wpn)|  q-(x')  qn(xn) |  | dx1 |  |dxn|  . 

1*1,  n 

For  n  even  there  will  be  one  more  term  and  the  designated  integral  is 


bounded  by 

(Hd))(n'2)/2  n  |  |  q1 1  |  f  |J  (wpQ)  | 
1*1, 2, n 


||H(wp1)||  q'(x')  q2(x2)  qn(xn) 
| dx1 |  | dx2 |  | dxn | 


But 

J  |  H(u>p1)  |  |  q2  (x2 )  |  |  dx2  |  <  c(  |  log  u>|  +  1)  ||q2|l 

and  for  w  > 


for  any  R.  Thus  we  have  for  n  even  or  odd,  the  bound  given  by  (3.1). 

To  complete  the  finding  of  the  bound  we  can  set  x  -  0.  Then  Pq  "  I*'  + 
x  | .  From  the  reglarity  properties  of  f  for  finite  argument  and  its 
behavior  at  «  we  can  estimate  the  remaining  integral  by  | |q^| |  | |qn| |  as  in 
the  previous  cases.  This  completes  the  proof  of  Lemma  3.2. 

This  leads  us  immediately  to 
Lemma  3.3.  Let 

*  w  n-1  n“*  n  i  i  i 

q  -  J  du  w  J (wpn)  (  n  H(wp.)  +  n  H(-wp.)  II  q  (x  )  |dx  |  . 
n  0  U  0  1  0  1  i-1 

*  n  i 

Then  |q  |  <  c  II  |  |q  |  |  for  n  >  5. 
n  i-l 

Proof:  From  the  decay  of  J,  H  as  w  -»  ®  we  see  that  q^  satisfies  the 

desired  bound  by  using  Lemma  3.2. 

n  i 

Lemma  3.4.  |q  |  Sc  II  ||q  ||  for  n  a  5.  Here  q  is  defined  in 
n  i-l 

(1.6). 

★  i  i  i 

From  (1.6)  and  (1.7)  it  is  easily  seen  that  q^  -  qn  if  q  (x  )  -  qix 
and  we  integrate  over  6  to  obtain  the  Bessel  function  as  we  did  in  deriving 
q^.  Hence  Lemma  3.4  follows  from  Lemma  3.2. 


-•  bounds  for  n  -  3,4. 

To  handle  the  misting  cas«t  n  •  3,4  we  must  use  s  less  crude  esci.i .,»• 
and  take  Into  account  the  oscillatory  character  of  the  Bessel  and  Hai.t'e. 
function.  However  it  is  such  easier  to  carry  this  out  for  n  -  3,4  chan  : 
n  -  2  because  the  integrand  is  dying  out,  if  slowly,  u  a» 

*  n  i 

T  4.1.  For  n  -  3,4,  |q^|  is  bounded  by  II  ||q  ||  and  jq^,  it 

bounded  | | q | | n. 

Proof:  The  integrand  of  defined  in  Leona  3.3  dies  out  respect ivriv 

-1/2  *  1 

like  u>  '  ,  u>  for  w  real  and  p^  >6  >0  say.  But  for  w  in  the  complt-x 

plane  with  a  suitable  argument  +  0,  q*  dies  out  exponentially.  We  split  thr 

Bessel  function  J(wp^)  as  +  H^^(-wPq)  for  real  and  note  that 

-  i  H^(u>p^)f  see  section  2.  Next  we  note  that  H^^(*p.<  - 

- 1/2 

exp  i  (u>p^  -  w/4)  (wp^)  for  large  |u>p^|,  |  arg  w|  <  6  say.  Hence  tor 

*  0  the  integrand  of  q*  behaves  for  large  |u>|  ,  up  to  a  phase  factor  aru:  t 

.'onscant,  like 
-1/2  -1/2 

~  (p0plp2)  (exp  iw(^o*pl*p2*  *  exp  lw( ‘po+pl+p2)  *  eXp  lw(p0'pi*"/ 

+  exp  iu)(-pQ-p1-p2)) 

for  n  -  3  and  a  similar  expression  appears  for  n  -  4.  But  the  factor  i  a 

-1,  x - 1/2 

( P0plp2p3^ 

First  we  note  that  the  integral  over  converges  in  the  neighborhood  a: 
j  *  0  is  it  did  for  n  i  5.  Without  loss  of  generality  we  consider  one 

the  exponential  terms  and  a  path  of  integration  in  the  w  plane  going  ti 
oela  as  arg  w  ■  a  to  infinity.  There  is  no  contribution  from  an  arc  at  •-  *: 
we  shift  to  this  path.  The  value  of  a  is  chosen  small  but  so  that 
exponential  is  decaying. 

The  easiest  case  is  the  first  term  for  n  -  4.  We  have  the  exponent 
pO*pl*p2*p3  >  0  80  a  >  S#t  w(p0^pl+p2+p3)  "  0n  inteBraCin&  wi:; 


respect  to  Q  wo  obtain  cip^p^p ^p^) 
:erms  wo  find 


Thus  doing  cho  iam  with  the  other 


4  4 

|q*|  Sen  | | ql | |  ♦  c  /  <P0P1P2P3)  1/2  ®  Nl(xl>l  ldxl| 

With  x  -  0  and  bounding  ^oplp2p3^  by  lxl+x^l  |x*-xz| 

Ux^-x3|  *  4-  |x3-x4|  S  wo  estimate  by  intogratlng  first  w.r.t  x3  in  the 

3  2 

usual  way.  This  glvos  tho  factor  | |q  ||.  Repeat  w.r.t.  x  and  again  w  r  r 
x*4  and  finally  w.r.t.  x*  Thus  wo  find 

*  4  i 

MV'-  5  c  n  m  u 

i 

★  co 

To  estimate  q^  in  L  wo  use  tho  same  idea  but  wo  obtain 


!q*|<cir|  | q 1 1  l+f(fi0/>lfi2)  'l//2(c(fiQ+fil+p2)'l//2+cl  -fiQ^pl*p2l  'l/2)  II  ql(xl)  |  dx  1  | 

•1/2  -1/2 

We  estimate  tho  first  term  using  (PqP^P 2)  (p0+pl+p 2^  bounded 

tj v  |x^+x3|  *(|x^-x^|  l*,x2.«3,  and  obtain,  as  for  n  -  4,  an  estimate  for 

3  t 

the  first  tera  of  tho  fora  c  IJ  /|q  //. 

1-1 

Tho  second  term  is  more  difficult  because  ‘^q*p1*p 2  vanishes  on  an 

el  1  ipse  if  wo  fix  x*  and  x3  this  follows  from  the  fact  that  p^+p ^  ls  -he 

2  13 

.uni  of  the  distances  of  x  to  x  and  x  .  First  we  use  six  dimensional  polar 

coordinates  to  estlaate  the  contribution  from  the  neighborhood  of  ^  ^  - 

2  2  2  2 

.  ,  -  0.  Thus  with  Pq  ♦  p^  +  p^  •  R  the  contribution  from  R  <  Rq  to  the 

unestimated  Integral  in  (4.1)  can  be  estimated  by 

II  ||q*||  /  |  -  cos  4  ♦  sin  ^  sin  ( $  ♦  ?)  |  d0  dR 

R<R0 

where  p^  -  R  cos  4,  P\  ~  R  sin  ♦  cos  ^2  "  R  sin  ^  s*n  **  T^ie  inteSraI 

converges  and  we  have  a  bound  with  the  factor  Rq . 

To  estimate  for  R  >  Rq  we  first  integrate  over  x"  in  the  region 

| -Pq+  p^+  p^\  >  6  and  obtain  using  bounds  n  q*  the  term 

C| |q2| |«‘1/2  /  Pol/2(qL(xl)|  |q3(x3)|  | dx1 |  |dx3|  . 
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since 

1  J  (plp2)  V2  q2(x2)  |dx21  <  f  (pl1+  P2X)  lq2(x2)l  dx2> 

<  c| | q2 | |  . 

Estimating  the  remaining  integrals  we  obtain 

C  S'1/2  *1 | q1 | |  . 

Over  the  neighborhood  of  the  ellipse  -p^+p^+p 2  *  0  i.e.  I ‘^o+^l+^2 ^  5  * 

we  use  L*  estimates  for  q*  and  estimate  the  area  by  the  width  6  times  the 

2 

length  of  the  ellipse  L.  Thus  we  obtain  from  the  integration  w.r.t.  x  , 

C  Rq1! |q2| I  f  ,01/2  S1/2  L  Iq1!  |q3|  Idx1!  |dx3|  , 

13  12 

where  L  is  the  length  of  the  ellipse  with  foci  at  x  ,x  and  with  |x  -x  |  + 

| X2 - X3 1  -  |x*+x^|.  Thus  L  is  bounded  by  (x^+x^l  -  p q.  Choosing  6  - 

l  3-1 

|x  +x  |  we  finally  have  the  bound  for  arbitrary  RQ  for  the  last  part: 

OR'1!  |q2|  |  ||q1||co||q3||-  <  cR*^  Mq1!!  . 

We  have  completed  the  proof  of  the  first  part  of  Lemma  4.1.  The  second 
part  follows  from  the  definition  of  in  (1.6)  and  (1.7). 


5 .  L 


n  £  3 . 


We  must  now  take  into  account  the  sign  of  q.  First  we  look  at  the 

multiple  kernel  that  occurs  in  q*.  It  is 

*  n-1  n-1 

K  -  X  <*>  j (wpft)  (  n  H(w p.)  +  n  H(-wp.)  )  dw  . 

0  U  l-l  1  i-1  1 

For  n  >  3  the  integral  converges  improperly.  We  claim 

I-gnmu.  S  1  K(  - 1)"  i.  0. 

Proof:  We  have  H(wp)  -  -c+  ^wt(t 2-p2)  *^2  dt  where  c+  >  0.  Also 

iJ  +  N  —  H  so  that  J  —  c+  S x>p  s^n  dt.  Thus  we  may  write 

,  «  n-1  -iwt.  n-1  iwt.  n  o  o  i  /o 

K  -  (-l)n'A  J  J  w  sin  wt- (  ne  1  +  II  e  l)n  (tf-pf)  '  dt.  dw. 

0  t,>p,  u  1  1  0  1  1  1 


0  t£>Pl  1  1  0  *  * 

The  integral  with  respect  to  u>  would  be 

*  t£)  iu>(-t0-S  t£)  iw(tQ+Z  t.)  ia)(-t0+Z  t£) 

ij(e  -e  +  e  -  e  )w  du> 


Re  iwe^01*7  dw  -  Re  f  e^^*7  du>  -  S' (a). 
0  da  0 


Hence 


&  -  (-Dn  C+  J  [«'(t0-  2  tt)  -  5  ’  ( -  tQ-S  t£)  +  <’(t0+  2  t£> 

-  5 ' ( - tg+  2  t ^ ) ]  n  (t2  -  p2)*1^2  dt^  for  tQ  >  pQ,  t^  >  p^ 

-  (-l)nc+(«'(t0-  2  tp  -  fi’(-t0-  2  tt)  )  n  (t2-p2)'1/2  dti 

for  tQ  >  p0,  tt  >  pt  . 

There  is  no  contribution  from  the  second  term.  Integrating  by  parts  wi: 


respect  to  t^  we  find 


K  -  (.1)"  C.  J  '0<‘0-'0>'3/2  "  • 

ti>pi>0  1 

where  t^  -  I  t^  since  t^  <  0  is  inadmissible.  Hence  the  kernel  has  a 


that  depends  only  on  n. 
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Corollary  to  La««^  £  1  alternates  In  sign. 

If  now  vt  write  q*  -  q*  -  q*  whara  ♦  (-)  indicates  the  positive  (and 
negative)  part  we  have  both  q*  and  q*  lying  in  Q  if  q*  does.  Then 

*  f.Il  1/^,1 

q  -  J  K  II  q  (x  )  dx 
n  1 

and 

S  lq*l  |dx|  S  S  K  n  (q*  ♦  ql)  dx1 
"  1 

-  Z  /  K  n  q|  dx1  . 

where  I  means  over  all  possible  choices  of  sign.  Hence  we  have: 

t  mmmM  5.2.  q*  is  bounded  in  L*  by  2n  II  ||q*||  if  /  K  II  |q*|  dx*  is 
n  1 

bounded  in  lA  by  II  |  |  q*  |  I  if  q*  >  0. 

T  *itm  s .  3 ,  jKIIq*dx*is  bounded  in  L*  by  II  ||q*||  if  q*  >  0,  n  >  J 
Proof:  Here 

«  n  n 

K  -  /  u>  J (wpn)  (  n  H (wp  )  ♦  n  H <-wp.)  )  dw 
0  U  i-1  1  i-i 

lu  obtain  the  estimate  we  may  Integrate  J(u>Pq)  with  respect  to 
1  tr  w.r.t.  Pq  -  |2x-x*-xn|.  We  use  ujJ(wpQ)  |dx|  -  wJ(wpQ)p0  dp()  do 

-J(WP0)P0)#  dp0  d#0 

We  anticipate  that  the  integral  converges  very  weakly  so  we  consider 

llqjl  -  L.  *  li«  L<L)  . 

L-*« 

.her* 

•  n  1  1  i 

l  -  f  du,  /  K  n  q  (x  )  | dx  I  | dx |  . 

L  *  pl<L  1 

VL 

•  n  n  n  ,  -  ■ 

U  .  f  du  r  lJ(wL)  (  n  )  ♦  n  H(-w* .))  n  q  (x  )  |dx  I 

L-«  1  1  1  1  1 


4 


•>he  re 


L9 


Split  tha  integral  over  u>  into  the  part  from  L  *  to  Qq  and  frm  to  ®. 
Consider  the  second  part  first  and  integrate  by  parts  noting  that  LJ (u>L)  du> 
-  dJ(u>L).  Thus,  taking  only  one  product,  since  the  other  one  is  essentially 
the  same, 

Xi  CO 

f  l5(o>l)  n  H(u>p . )  dw  -  -  j(nnL)  n  mo.p  )  -  J  j(wL)  s  p.fi<wp  >  n  H(wp.)dw. 

%  %  i  1  J-i  1 

The  first  term  Is  bounded  for  OqP^  k,  L'^Oq1/2!!  (|log  Cl0|  +  |log  p x\ ) . 


If  we  integrsts  over  II  q*  dx*  and  use  J  |  log  p^|  ql  dx  1  <  | |ql| |  and  over 

n  q*  dx*  where  p^  >  const.  and  H(OqP^)  is  bounded  we  find  that  the 

- 1/2  i 

integrated  term  dies  out  at  least  like  L  /  n  | |q  | | .  Thus  we  are  reduced 
to  bounding  for  large  L, 

«0 

K  -  /  J (u>L)  2  p.  fi(wp. )  n  H(wp.)  dw  . 

wo  n0  i  1  1  j*i  1 

.»pp lied  to  n  q  *  dx*  . 

Collapsing  the  integral  using  Lemma  3 . 1  we  find  a  typical  term  (i  -  1> 

’•ounded  by 

n  • 

n  |  1  q 1 1  |  /  |  J(u>L)  |  p.|ft(up  )  q* (x* )  |dx*|  fl(w)(n‘2)/2  du, 

1-2  % 


i  ^  ,,.t, 


/  pl|fl(wpl)|ql(x1)  |  dx1 1  <  Hq1!!*’1  ♦  f  (upl)'L/‘i  px  q1  |  dx 


1.  ,  -1 


.-V2  „  l  1 


.t<L  wpL>k 

^  . .  1, ,  ,  -1  A  , 1/2  -1/2 . 

<  I |q  1 1  (w  ♦  L  '  w  ). 

nence  for  large  L  the  bound  for  the  integral  from  0^  to  *  is 

n  )  | q 1 1  |  /  (U.L)'1/2  (w  1  L1/2  u»'1/2)  H(ui)  2  dw 

0„ 


<  n  | | q1 | |  (L'l/2  log  0Q2  0q1/2  *  con,t  >  for  n  s  3 


Sen  | |q  || 


for  n  z  3 


.u:e  that  this  term  does  not  go  to  zero  as  L  -•  • 


The  same  estimate  holds  for  the  remaining  term. 

We  are  left  with  the  integral  from  L€  to  Oq.  This  is  the  trickiest 
and  reflects  the  poor  decay  properties  of  the  time  dependent  problem,  see 
section  1. 

For  this  integral  we  can  consider  L  sufficiently  large  so  that  J(u>L)  - 
(wL)  ‘1/2  cos(wL  -  J)  .  Setting  <*>L  -  0,  the  kernel ,  which  we  shall  show  tends 
to  0  is 

°0L 

f  n'1/2  sin(o  -  J)  (n  H(n<7i)  +  n  hc-o^))  dn  , 

L1’  * 

where  a ^  -  p^/L  <  1 . 

If  >  6  for  at  least  two  values  of  i  the  limit  of  this  integral  as  L 
-  jo  would  be  zero  because  of  the  decay  of  the  Hankel  functions. 

We  consider  the  integral  in  the  complex  0  plane  on  a  path  chosen  so 
that  the  integral  Is  decaying  exponentially.  Thus  a  typical  term  is  the 
integral  of  ft  ^2exp  i(0-*/4)  II  H(Oo^)  on  |0|  -  L*  ,  0  -  e^a|Q|  and  |W|  - 
4i  L.  The  integral  on  the  ray  for  a  *  0  decays  exponentially  and  is  bounded 
bv 

/  |0| ' 1/2  e'a|Q|  n  |H(|0|<7.  ela)|  d|Q| 

L1  *  * 

<  c  f  |Q|  ‘ 1/2  e‘a|fl|  n  (|log  OaJ+1)  d|0| 

<  c  L'(l'<)/2  e‘aLl  a" L( | log  L* |  +  log  p£)  . 

•  iitegrated  with  respect  to  n  q^(xb  dx^  we  find  that  this  term  -»  0  as  L  -  u 
tor  | | q 1 1 |  bounded. 

The  integral  on  |ft|  -  L*  la  bounded  by 

c  l1/2  f  e * 1  *in  9  n  H (to.  e1^ )  e19  At  , 

0 

.  ith  l  -  d  *  or  for  small  a  by 


f 


n 


q ( X )  +  Tq  -  qR  , 

where  Tq  is  a  convergent  sum  for  ||q||  sufficiently  small,  bounded  by 
2 

| | q | |  Hence  if  q^  €  Q  and  is  sufficiently  small  this  equation  can  be 

solved  for  q  uniquely. 


Appendix  1.  Derivatin  of  (1.4)  and  (1.4)*. 

To  derive  (1.4)*  we  use  the  integral  equatin  derived  from  (1.2)*  along 
with  Pg  *  0  by  using  the  2D  Riemann  function.  Again  c  is  a  given  constant 
(Al.l)  Ps  -  c  /  ((t-t')2  -  |x-x'|2)'1/2  («(t'-e.x') 

+  Ps(e,t' ,x' ))  q(x')  |dx' |dt' 

The  far  field  in  the  back  scattered  direction  -e  is  given  by  S(r,e)  and 
is  obtained  from  Pg  by  setting  t-|x|  -  r,  x  *  -e|x|,  t+|x|  ®  which  yields 

c  /  ( r  - 1'  -e*x' )  q(x')  (6(  t'-e«x')  *  P^e.t'.x'))  |dx'|  dt  for  S  (r,e)  or 

c  /  (r-a)  q(x')  (5(a-2e»x')  +  Ps(e ,a-e»x' ,x' )  da  |dx'| 

This  is  the  Abel  transform  of 

/  q(x# ) (6 (a-2e»xl )  *  P^fe.a, e*x' ,x' ) )  |dx'|  . 

So  we  may  regard  this  quantity  as  given  by  the  back  scattered  field.  The 
first  term  is  of  lower  order  in  q  and  is  also  a  Radon  transform  of  q.  Hence 
if  we  regard  the  Born  approximation  as  given  R.T.  q^  -  R.T.  q  +  f  q(x' ) 
P  (e,a-e*x' ,x' )  |dx'|.  Inverting  the  Radon  transform  we  obtain  (1.4)*.  To 
obtain  (1.4)  we  insert  (1.8)  and  integrate  with  respect  to  a.  Noting  that 

+jo  i  (jxj 

-  da  changes  sign  with  u>  for  real  but  is  independent  of  to  for  u  >  0 

-  00 

or  u)  <  0  we  easily  obtain  (1.4). 
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Appendix  2.  Convergence  of  an  Iterative  Procedure. 


q  +  Tq  -  qfi 

we  may  use  an  iterative  or  power  series  expansion  to  find  its  solution.  In 
[1]  we  have  attempted  the  iteration.  Here  we  show  the  iteration  converges 
for  IMgll  small.  The  method  is  to  find  at  the  iteration  from 

qN  “  qB  -  T  qN-l  • 

Thus  Wl  -  'TqN-l  +  TqN-2' 

Let  Q  -  max  | |qT| | ,  i.e.  the  max  norm  that  occurs  up  to  and  including 
W  JsN-1  J 

the  NC^  iteration.  Then  by  lemmas  2.2,  3.3,  4.1  and  5.2  if  we  define  q^  in 
the  obvious  way, 

""ns  •  vm"  s  ■  w  • 


Hence 


lqN  •  qN- 1  * 


We  also  have  with  q^^  -  q^ 


2  * ^  qN- 1 *qN- 2  I  • 


i  a  *  n_n 

lqnNM  S  c  Vl  ' 


so  that 


2.2 

iqN'i 5  iiqBn  r 


c2qn-i 

qn  +  • 

We  claim  if  -  ilq^M  sufficiently  small  the  iteration  converges. 

To  converge  we  need,  setting  cQ^  - 


<£  1  . 


(1-V‘ 


Claim.  If  <  l/16c  then  <  l/2c.  We  may  assume  c  >  10  say. 
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Proof:  By  the  inequality  for  we  have  if  <  l/2c, 

S  X.  *  £  Tf  *  - 1 - 

N  1  1'XN-1  16c 


4C(1  '  2c^ 


■^(i  + 


D  • 


2c  8  2(2c- 1) 

1  cl 

Hence  <  Jc  Provic*es  (2c  1)  <  4  or  <  ,  or  c  >  7/10. 


this  guarantees  that 


*N 


(1*V 


2  <  1.  Hence  the  iteration  converges, 
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